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Q^ ' Abstract 



Prediction and prevention of musculo-skeletal injuries is an important aspect of preventive 
health science. Using as an example a human knee joint, this paper proposes a new coupled- 
loading-rate hypothesis, which states that a generic cause of any musculo-skeletal injury is a 
Euclidean jolt, or ^^(S)— jolt, an impulsive loading that hits a joint in several coupled degrees- 
^^ ^ of-freedom simultaneously. Informally, it is a rate-of-change of joint acceleration in all 6-degrees- 

^^ . of-freedom simultaneously, times the corresponding portion of the body mass. In the case of a 

^\ ' human knee, this happens when most of the body mass is on one leg with a semi-flexed knee - 

and then, caused by some external shock, the knee suddenly 'jerks'; this can happen in running, 
skiing, sports games (e.g., soccer, rugby) and various crashes/impacts. To show this formally, 
based on the previously defined covariant force law and its application to traumatic brain in- 
pH , jury [Ivancevic 2008J . we formulate the coupled Newton-Euler dynamics of human joint motions 



and derive from it the corresponding coupled 5'-B(3)— jolt dynamics of the joint in case. The 
S'i5(3)— jolt is the main cause of two forms of discontinuous joint injury: (i) mild rotational discli- 
i-Q ■ nations and (ii) severe translational dislocations. Both the joint disclinations and dislocations, 

^^' as caused by the ^^(S)— jolt, are described using the Cosserat multipolar viscoelastic continuum 

joint model. 
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1 Introduction 



In this paper, we propose a new model of musculo-skeletal injury, using as an example the human 
knee joint. The knee joint comprises three articulations: (i) a tibio-femoral joint) between the medial 
and lateral condyles of the femur and tibia (see Figure [J), (ii) patelo- femoral joint between the 
femur and patella, and (iii) tibio-fibular joint between tibia and fibula. The knee is double condyloid 
joint, with a dominant flexion/extension. As a synovial joint, the knee has strong fibrous capsule 
that attaches superiorly to the femur and inferiorly to the articular margin of the tibia. Given its 
relatively poor bony fit, the knee relies on ligaments for much of its structural stability and integrity. 
[Whiting and Zernicke 1998 . 
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Figure 1: Schematic latero-frontal view of the left knee joint. Although designed to perform mainly 
flexion/extension (strictly in the sagittal plane) with some restricted medial/lateral rotation in the 
semi-flexed position, it is clear that the knee joint really has at least 6 degrees-of-freedom, including 
3 micro-translations. The injury actually occurs when some of these microscopic translations become 
macroscopic, which normally happens only after an external jolt. 

Knee joint injuries range from mild ligament or meniscus tearing to severe traumatic dislocations 
(see [Seroyer et al 2008] and references therein) that fall among the most severe form of ligament 
injury to the lower extremity, associated with a high rate of complications including amputation. 

Knee joint injuries are also frequent in sports, especially in ball games. For example, in a recently 
case-reported complex knee injury in a rugby league player [Shilhngton et al 2008| , resulting in com- 
bined rupture of the patellar tendon, anterior cruciate and medial collateral ligaments, with a medial 
meniscal tear, the video-analysis suggests two points during the tackle when there was the potential 
for injury: the first occurred when the player was in single leg stance whilst running and received 
impact to his upper body from three defenders; the second occurred when the player landed on the 
knee and sustained a valgus and hyper-flexion force under the weight of two defenders. The goal of 
treatment in this condition is restoration of both the extensor mechanism and knee stability. 

Also, the increased number of women participating in sports like soccer has been paralleled 
by a greater knee injury rate in women compared to men. In particular, menstrual cycle phase 
has been correlated with risk of noncontact anterior cruciate ligament injury in women (see, e.g. 
[Chaudhari et al 2007| )Pl Among these injuries, those occurring to the anterior cruciate ligament 
are commonly observed during sidestep cutting maneuvers [Sanna and O'Connor 2008J . In addition, 
general fatigue appears to correlate with injuries to the passive knee-joint structures during a soccer 



^According to Chaudhari et al 2007'. jumping and landing activities performed during different phases of the men- 
strual cycle lead to differences in foot strike knee flexion, as well as peak knee and hip loads, in women not taking an 
oral contraceptive but not in women taking an oral contraceptive. Women will experience greater normalized joint loads 
than men during these activities. 



game. It has been observed that a higher injury rate and more severe injuries occur towards the end of 
a soccer game or practice [ Hawkins et al 200T1 |Ostenberg and Roos 2000| , suggesting a fatigue effect 
on the neuromuscular system. Besides, women soccer players who sustain knee injuries have a high 
risk of going on to develop osteoarthritis at a young age [Lohmander et al 2004] . 

Besides, the knee is the body part most commonly injured as a consequence of collisions, falls, and 
overuse occurring from childhood sports. The number of sports-related injuries is increasing because 
of active participation of children in competitive sports [Slow et al 2008| . Children differ from adults 
in many areas, such as increased rate and ability of healing, higher strength of ligaments compared 
with growth plates, and continued growth. Growth around the knee can be affected if the growth 
plates are involved in injuries [Slow et al 2008J . 

Literature on knee mechanics has been reviewed in |Komistek et al 2005] , evaluating various tech- 
niques that had been used to determine in 'vivo loads' in the human knee joint. Two main techniques 
that had been used were 3-axial accelerometer-based telemetry - an experimental approach, and math- 
ematical modelling - a theoretical approach. Accelerometric analyzes had previously been used to 
determine the 'in vivo' loading of the human hip and more recently evaluated in the determination 
of in vivo knee loads. Mathematical modelling approaches can be categorized in two ways: (i) those 
that use optimization techniques to solve an indeterminate system, and (ii) those that utilize a reduc- 
tion method that minimizes the number of unknowns, keeping the system solvable as the number of 
equations of motion are equal to the number of unknown quantities (for more technical details, see 
[Komiste k et al 2005] and references therein). 

The use of a force-controlled dynamic knee simulator to quantify the mechanical performance of 
total knee replacement designs during functional activity was pioneered by [DesJardins et al 2000] . 
in which dynamic total knee replacem,ent (TKR) study utilized a 6-degree-of-freedom force-controlled 
knee simulator to quantify the effect of TKR design alone on TKR mechanics during a simulated 
walking cycle. Simultaneous prediction of implant kinematics and contact mechanics has been demon- 
strated using explicit finite element (FE) models of the Instron/Stanmore Knee Joint Simulator, In- 
stron, Canton, MA |Godest et al 2002], IHalloran et al 2005a| IHalloran et al 2005bj In these models, 
kinematic verification was performed by comparing experimental and model-predicted motion for a 
single implant. Both models were found to produce similar kinematic simulation results. Estimated 
contact pressure distributions were also closely correlated, as long as significant edge-loading con- 
ditions were not present. Recently, an adaptive FE method for pre-clinical wear testing of TKR 
components was developed in [Knight et al 2007] , capable of simulating wear of a polyethylene tib- 
ial insert and to compare predicted kinematics, weight loss due to wear, and wear depth contours 
to results from a force-controlled experimental knee simulator. The displacement-controlled inputs, 
by accurately matching the experimental tibio-femoral motion, provided an evaluation of the simple 
wear theory. The force-controlled inputs provided an evaluation of the overall numerical method by 
simultaneously predicting both kinematics and wear. Proposed international standards for TKR wear 
simulation have been drafted (see [ISO Standard 14243-2, 2000] ), yet their methods continue to be 
debated |Laz et al, 2005| . The 'gold standard' to which all TKR wear testing methodologies should 
be compared is measured in vivo TKR performance in patients. The study of [DesJardins et al 2007J 
compared patient TKR kinematics from fluoroscopic analysis and simulator TKR kinematics from 
force-controlled wear testing to quantify similarities in clinical ranges of motion and contact bearing 
kinematics and to evaluate the proposed ISO force-controlled Stanmore wear testing methodology. 

For human movement purposes, we can say that the safe knee motions (flexion/extension with 
some medial/lateral rotation in the flexed position) are governed by standard Euler's rotational dy- 
namics coupled to Newton's micro-translational dynamics. On the other hand, the unsafe knee events, 
in this paper we will show that the main cause of knee injuries are the knee SE(3)-jolts, the sharp and 
sudden, "delta"- (forces -f- torques) combined. These knee SE(3)-jolts do not belong to the standard 
Newton-Euler dynamics. The only way to monitor them would be to measure "in vivo" the rate of 
the combined (forces -I- torques)- rise in the knee joint (see Figure [T]). 

This paper proposes a new hypothesis for generic musculo-skeletal injury, called: 

Coupled loading rate hypothesis: 

The main cause of knee injury is Euclidean jolt, or SE{3)—jolt, an impulsive loading that hits the 
knee joint in several coupled degrees- of-freedom (DOF) simultaneously. The same hypothesis applies 



to any other major human joint: all acute musculo- skeletal injuries are caused by some form of Eu- 
clidean jolt. 

In the realm of musculo-skeletal injury, a Euclidean jolt represents a 6-degree-of-freedom 'jerk' 
(rate-of-change of acceleration) times most-of-tlic-body mass. In other words, a Euclidean jolt is a 
time derivative of the Euclidean force (three-dimensional force + three-dimensional torque). In the 
case of a knee, it happens when all the body mass is on one leg with a semi-flexed knee - and then, 
caused by some external shock, the knee 'jerks'; this often happens in running, skiing, sports games 
(e.g., soccer, rugby) and various crashes/impacts. The reason why this happens in this particular 
scenario is the following: in a semi-flexed position, the knee joint has all 6 DOFp if a sudden jerk 
with the full body mass hits the knee in this position, it will happen in all 6 DOF simultaneously, 
tearing apparat the soft tissue. If the resulting jolt is of high intensity then even the hard tissue will 
be damaged. 

To demonstrate this formally, based on the previously deflned covariant force law, we formulate 
the coupled Newton-Euler dynamics of the knee motions and derive from it the corresponding coupled 
S* £'(3)— jolt dynamics. The effect of 5'i?(3)— jolt can be seen in two forms of discontinuous knee injury: 
(i) mild rotational disclinations and (ii) severe translational dislocations. Both the knee disclinations 
and dislocations, as caused by the ^-^(B)— jolt, are described using the Cosserat multipolar viscoelastic 
continuum model. 

While we can intuitively visualize the knee SE(3)-jolt, for the purpose of whole human musculo- 
skeletal dynamics simulation, to avoid deterministic chaos caused by nonlinear coupling, we use the 
necessary simplified, decoupled approach (neglecting the 3D torque matrix and its coupling to the 3D 
force vector). In this decoupled framework of reduced complexity, we define: 

The cause of knee dislocations is a linear 3D-jolt vector, the time rate-of-change of a 3D-force 
vector (linear jolt = mass x linear jerk). The cause of knee disclinations is an angular 3-axial jolt, 
the time rate-of-change of a 3-axial torque (angular jolt = inertia moment x angular jerk). 

This decoupled framework has been implemented in the Human Biodynamics Engine [Ivancevic 2005J , 
a high-resolution neuro-musculo-skeletal dynamics simulator (with 270 DOFs, the same number of 
equivalent muscular actuators and two-level neural reflex control), developed by the present author 
at Defence Science and Technology Organization, Australia. This kinematically validated human 
motion simulator has been described in a series of papers and books fivanccvic and S noswell 2001[ 
[Ivancevic and Beagley 20031 Hvancevic 2002( llvancevic 2004) |lvancevic~a nd Beagley 2005 , 
[Ivancevic and Ivancevi c 2006a, Iva ncevic and Ivancev ic 2006b, Ivancevic and Ivancevic 2006c], 
[Ivancevic and Ivancevic 2007d, Ivancevic and Ivancevic 2007c, Iva ncevic 200 6, Ivancevic and Ivancevic 2007a| 
llvancevic and Iv ancevic 2006. Ivancevic and Ivancevic 2007 b. Iva ncevic and Ivancevic 20(l8] . 

2 The S'£'(3)— jolt: the main cause of human joint injury 

In the language of modern biodynamics [Ivancevic 2004[ llvancevic and Ivancevic 2006aj , the general 
knee motion is governed by the Euclidean SE(3)-group of 3D motions. Within the knee SE(3)-group 
we have both SE(3)-kinematics (consisting of the knee SE(3)-velocity and its two time derivatives: 
SE(3)-acceleration and SE(3)-jerk) and the knee SE(3)-dynamics (consisting of SE(3)-momentum 
and its two time derivatives: SE(3)-force and SE(3)-jolt), which is the knee kinematics x the knee 
mass-inertia distribution. 

Informally, the knee SE(3)-jolv^ is a sharp and sudden change in the SE(3)-force acting on the 



^The frequently used human leg model by IBrand et al 1994] was comprised of 47 muscles and each joint was repre- 
sented by three interactive forces and three interactive torques jK omist ek et al 200 5 - thus clearly showing a generic 
6 DOF function. In the case of the knee joint, the macroscopic movement of patella during flexion/extension clearly 
shows the existence of micro- translations in this joint (see also explanation in the caption of Figure 1). As this paper 
tries to show generic musculo-skeletal injury patterns, note that similar 6 DOF function is even more prominent in the 
case of human shoulder. 

^The mechanical SE(3)-jolt concept is based on the mathematical concept of higher— order tangency (rig- 
orously defined in terms of jet bundles of the head's configuration manifold) [Ivancevic and Ivancevic 2006cl 
llvancevic and Ivancevic 2007e] . as follows: When something hits the human head, or the head hits some external 
body, we have a collision. This is naturally described by the SE(3)— momentum, which is a nonlinear coupling of 3 linear 
Newtonian momenta with 3 angular Eulerian momenta. The tangent to the SE(3)— momentum, defined by the (absolute) 
time derivative, is the SE(3)— force. The second-order tangency is given by the SE(3)— jolt, which is the tangent to the 



mass-inertia distribution within the knee joint. That is, a 'delta'-change in a 3D force-vector coupled 
to a 3D torque- vector, hitting the knee. In other words, the knee SE(3)-jolt is a sudden, sharp and 
discontinues shock in aU 6 coupled dimensions of the knee joint, within the three Cartesian (x, y, z)- 
translations and the three corresponding Euler angles around the Cartesian axes: roll, pitch and yaw 
[Ivancevic and Beagley 2003| . If the SE(3)~jolt produces a mild shock to the knee (internal sudden 
loss of stability), it causes mild, soft-tissue knee injury. If the SE(3)-jolt produces a hard shock 
(external hit by a massive body) to the knee, it causes severe, hard-tissue knee injury, with the total 
loss of knee movement. 

The knee SE(3)-jolt is rigorously defined in terms of differential geometry 
[Ivancevic and Ivancevic 2006cllTvancevic and Ivancevic 2007e| . Briefly, it is the absolute time-derivative 
of the covariant force 1-form (or, co- vector field) applied to the knee. With this respect, recall that the 
fundamental law of biomechanics - the so-called covariant force law [Ivancevic and Ivancevic 2006bi 
llvancevic and Ivancevic 2006c[ Hvancevic and Ivancevic 2007e] . states: 

Force co-vector field = Mass distribution x Acceleration vector-field, 

which is formally written (using the Einstein summation convention, with indices labelling the three 
local Cartesian translations and the corresponding three local Euler angles): 



F„ 



rriu^a 



{fi, V — 1, ..., 6 = 3 Cartesian + 3 Euler) 



where F^ denotes the 6 covariant components of the knee SE(3)-force co-vector field, m^j^ represents 
the 6x6 covariant components of the inertia-metric tensor of the total mass moving in the knee joint, 
while a^ corresponds to the 6 contravariant components of the knee SE(3)-acceleration vector- field. 

Now, the covariant (absolute, Bianchi) time-derivative |j(-) of the covariant SE(3)-force F^ defines 
the corresponding knee SE(3)-jolt co-vector field: 



D 

~dt 



(F^) = '^M-§(«') = '"A"- («" + ^W^^) 



(1) 



where |j(a'^) denotes the 6 contravariant components of the knee SE(3)-jerk vector-field and overdot 
(') denotes the time derivative. FJ^^ are the Christoffel's symbols of the Levi-Civita connection for the 
SE(3)-group, which are zero in case of pure Cartesian translations and nonzero in case of rotations 
as well as in the full-coupling of translations and rotations. 

In the following, we elaborate on the knee SE(3)-jolt concept (using vector and tensor methods) 
and its biophysical consequences in the form of the knee dislocations and disclinations. 

2.1 S'i?(3)— group of local joint motions 

Briefly, the S'_E(3)— group of knee motions is defined as a semidirect (noncommutative) product of 3D 
knee rotations and 3D knee micro-translations, 

SEii) :=5'0(3)[>R^ 



Its most important subgroups are the following: 



Subgroup 



Definition 



5*0(3), group of rotations 
in 3D (a spherical joint) 



Set of all proper orthogonal 
3 X 3 — rotational matrices 



SE{2), special Euclidean group 
in 2D (all planar motions) 



Set of all 3 X 3 — matrices: 

cos 9 sin 6 r^ 

— sin 9 cos 9 r„ 











'y 
1 



SO{2), group of rotations in 2D 

subgroup of S'£'(2) -group 

(a revolute joint) 



f^ group of translations in 3D 
(all spatial displacements) 



Set of all proper orthogonal 
2 X 2 — rotational matrices 
included in SE{2) — group 



Euclidean 3D vector space 



SE(3)-force, also defined by the time derivative. 



In other words, the gauge S* £'(3)— group of knee Euchdean micro-motions contains matrices of the 



form 



R p 



where p is knee 3D micro-translation vector and R is knee 3D rotation matrix, given 



by the product R = R^p-R^-Re of the three Eulerian knee rotations, roll = R^p, pitch — R^, yaw — Rg, 
performed respectively about the a— axis by an angle ip, about the y— axis by an ang le ip, and about 
the z— axis by an angle 9 (see [Ivancevic 20041 [Park and Chung 2005[ llvancevic 2006| ). 



Ri^ 



1 

COS tfi — sin ip 
sin if cos ip 



Rd 



cos ip sin tjj 

1 

— sin ip cos tp 



Re 



cosf 

sin^ 





- sin e 

cos 61 

1 



Therefore, natural knee S'_E(3)— dynamics is given by the coupling of Newtonian (translational) 
and Eulerian (rotational) equations of the knee motion. 

2.2 Local joint 5" £"(3)— dynamics 

To support our locally-coupled loading-rate hypothesis, we formulate the coupled Newton-Euler 
dynamics of the knee motions within the S'i5( 3)— group. The forced Newton-Euler equations read in 
vector (boldface) form 



Newton 
Euler 

where x denotes the vector cross product. 



p = Mv = F + p X a;, 

•;T=ItJ = T + 7rxa; + pxv, 



(2) 



M = Mij — diag{mi, m2, rris} and 



diag{Ii, 12,13}, (i,j = 1,2,3) 



are the total moving segment's (diagonal) mass and inertia matricesp defining the total moving 
segment mass-inertia distribution, with principal inertia moments given in Cartesian coordinates 
{x, y, z) by volume integrals 



h = 



p{z'^ + y^)dxdydz, I2 = 



p{x'^ + z'^)dxdydz, I^ = 



p(x^ + y^)dxdydz, 



dependent on the knee density p — p{x, y, z). 



h;i,U2,W3 



and 



[uJi,UJ2,LlI3\ 



(where [ ]* denotes the vector transpose) are linear and angular knee-velocity vectors (that is, column 
vectors) , 

F = F, = [Fi, F2, F3] and T = T, = [T,,T2, T3] 

are gravitational and other external force and torque co-vectors (that is, row vectors) acting on the 
knee, 



p = pi = Mv ^[pi,p2,p3] = [mivi,m2V2,m2V2] 

TT = TTi = Iw = [tTi, 7r2, TTg] = [IiUJi, l2(^2, h^^s] 

are linear and angular knee-momentum co-vectors. 

In tensor form, the forced Newton-Euler equations ^ read 



and 



Pi 



M„?)J = F,, 



i+£^fcPJ^^ 



(z,i,/c = 1,2,3) 



*In reality, mass and inertia matrices (M,I) are not diagonal but rather full 3x3 positive-definite symmetric 
matrices with coupled mass- and inertia— products. Even more realistic, fully-coupled mass— inertial properties of a 
moving segment are defined by the single non-diagonal 6x6 positive-definite symmetric mass-inertia matrix MsE(3), 
the so-called material metric tensor of the 5-5(3)— group, which has all nonzero mass— inertia coupling products. However, 
for simplicity, in this paper we shall consider only the simple case of two separate diagonal 3x3 matrices (M, I). 



where the permutation symbol e^,^ is defined as 



k 

+1 if (i,j,A:) is (1,2, 3), (3, 1,2) or (2, 3,1), 

ik 



eik = <-l if (i, J, k) is (3, 2, 1), (1, 3, 2) or (2, 1,3), 



otherwise: i — j or j = k ot k — i. 

In scalar form, the forced Newton-Euler equations ^ expand as 

Pi = Fi — m3V3UJ2 + m2V2UJz 
Newton : l^ p^ ^ F2 + mswawi — miwicjs , (3) 

p^ = F3 — m2V2UJl + miViLU2 

TTj = Ti + (m2 - ni3)v2V3 + {I2 - 13)^2^3 
Euler : <( tt^ = T2 + (777,3 - 777i)7;i7;3 + (13 - Ii)ujiUJ3 , 

TT,, = T3 + (777i - m2)viV2 + {h - -^2)^1^2 

showing the moving segment's mass and inertia couplings. 

Equations Q^® can be derived from the translational + rotational kinetic energy of the moving 
segmenllfl 

Ek - -v*Mv + -lo'Ilo, (4) 

or, in tensor form 

For this we use the Kirchhoff-Lagrangian equations (see, e.g., [Lamb 19321 [Leonard 1997] , or the 
original work of Kirchhoff in German) 

^^d^Ek = d^Ek XLu + F, (5) 

-9^£;fc = d^Ek XUJ + d^Ek X V + T, 
at 

where dvEk = ^§^, d^Ek — -§^', in tensor form these equations read 

j^d^.E = el^{d,.E)uj^ + F,, 

j^d^.E = el,(d^,E)J^ + el,{d^,,E)v^^T.,. 
Using ([3])-([n]), linear and angular knee-momentum co- vectors are defined as 







p = (9v£^fc, 


TT = du^Ek, 


or, in tensor form 












Pi = d^iE, 


-Ki = d^iE, 


with their corresponding time derivatives, in vector form 






at at 


at at 


or, in tensor form 




Pt = -irPi = -jrOyiE, 
at at 


at at 


or, in scalar form 








P = [Pi 


,P2,P3' 


1 = [777lt)l,7772W2,"^3i'3], 


5T = [^l>2,Tr3] = [^1 



While healthy knee 5i?( 3)— dynamics is given by the coupled Newton-Euler micro-dynamics, the 
knee injury is actually caused by the sharp and discontinuous change in this natural SE{3) micro- 
dynamics, in the form of the S'-E(3)— jolt, causing discontinuous knee deformations, both translational 
dislocations and rotational disclinations. 

""In a fully— coupled Newton-Euler knee dynamics, instead of equation JJjl we would have moving segment's kinetic 
energy defined by the inner product: 



2.3 Joint injury dynamics: the S'£'(3)— jolt 

The 5-^(3)— jolt, the actual cause of the knee injury (in the form of the plastic deformations), is 
defined as a coupled Newton+Euler jolt; in (co)vector form the S* £'(3)— jolt readqj 

SE(3) - "oh ■ I Newton joh j F==p-pxw-pxw, 

\ Euler jolt : T — n —tt xw — TrxtJ — pxv — pxv, 

where the linear and angular jolt co- vectors are 

F = Mv = [Fi,F2,F3], T = Iu;= [Ti,T2,T3], 

where 

V = [iil,W2,W3]*, W = [aJi,(:02,W3]*, 

are linear and angular jerk vectors. 

In tensor form, the S'-E(3)— jolt readqj 

F^ = P^ - slkPo"^^ - e^fcPJ^^ ih 3, k ^ 1,2, 3) 

T^ = ^^ -^Ik^j^'' -^Ik^j^'' -^kPiv"" -e{^pjv'', 

in which the linear and angular jolt covcctors are defined as 

F = F, = Mv = Khj i^ = [Fi , F2 , Fs] , 
t = t, = Ii3 = hju^ = [Ti, t2, Ts], 

where v = ii', and w = w' are linear and angular jerk vectors. 
In scalar form, the 5'i?(3)— jolt expands as 

Fi=pi~ m2UJ3V2 + ms (^2^3 + V3UJ2) - m2V2ti>3, 
Newton jolt : ^ F2 — P2 + "t-i'^s'^i — ra^uoiv^ — m^vzioi + niiviuj^, 

F3 ^PS - miL02Vl + TO2'^l'^2 - V2LO1 - miViU!2, 

Ti = TTi - (m2 - ma) (W3W2 + W2i'3) - {h - I3) {UJ3UJ2 + ^2^3) , 

Euler jolt : <( t2 = tt2 + {mi ~ ms) {V3V1 + V1V3) + {h - I3) {uj^U!! + ujiujs) , 

T3 =TT3 - {mi - TO2) (w2t'i + wii'2) - (-^1 - h) {i-^20Ji + ^1^2) ■ 

We remark here that the linear and angular momenta (p, tt), forces (F, T) and jolts (F, T) are co- 
vectors (row vectors), while the linear and angular velocities (v, w), accelerations (v, Co) and jerks (v, w) 
are vectors (column vectors). This bio-physically means that the 'jerk' vector should not be confused 
with the 'jolt' co-vector. For example, the 'jerk' means shaking the head's own mass-inertia matrices 
(mainly in the atlanto-occipital and atlanto-axial joints), while the 'jolt'means actually hitting the 
head with some external mass-inertia matrices included in the 'hitting' SE(3)-jolt, or hitting some 
external static/massive body with the head (e.g., the ground - gravitational effect, or the wall - 
inertial effect). Consequently, the mass- less 'jerk' vector represents a (translational-l- rotational) non- 
collision effect that can cause only soft knee injuries, while the inertial 'jolt' co-vector represents a 
(translational-t-rotational) collision effect that can cause hard knee injuries. 

2.4 Joint disclinations and dislocations caused by the S" £'(3)— jolt 

For mild knee injury (caused by internal loss of stability), the best injury predictor is considered to 
be the product of localized knee strain and strain rate, which is the standard isotropic viscoelastic 
continuum concept (see, e.g. [Halloran et al 2005a| ). To improve this standard concept, in this sub- 
section, we consider the knee joint as a 3D anisotropic multipolar Cosserat viscoelastic continuum 



^Note that the derivative of the cross-product of two vectors follows the standard calculus product-rule: -jr (u X v) = 
u X V -I- u X V. 

In this paragraph the overdots actually denote the absolute Bianchi (covariant) time-derivative jlj, so that the jolts 
retain the proper covector character, which would be lost if ordinary time derivatives are used. However, for the sake 
of simplicity and wider readability, we stick to the same overdot notation. 



[Cosserat and Cosserat 18981 ICosserat and Cosserat 19091 [Eringen 2002] , exhibiting coupled-stress- 
strain elastic properties. This non-standard continuum model is suitable for analyzing plastic (ir- 
reversible) deformations and fracture mechanics [Bilby and Eshelby 1968| in multi-layered materials 
with microstructure (in which slips and bending of layers introduces additional degrees of freedom, 
non-existent in the standard continuum models; see ^Mindlin 19651 ILakes 1985J for physical charac- 
teristics and Yang and Lakes 1981HYang and Lakes 1982| , 
[Park and La 



tes 19861 for biomechanical applications). 



The £'£'(3)— jolt (F,T) causes two types of localized knee discontinuous deformations: 

1. The Newton jolt F can cause severe micro-translational dislocations, or discontinuities in the 
Cosserat translations; 

2. The Euler jolt T can cause mild micro-rotational disclinations, or discontinuities in the Cosserat 
rotations. 

For general treatment on dislocations and disclinations related to asymmetric discontinuous defor- 
mations in multipolar materials, see, e.g., |Jian and Xiao-ling 1995HYang et al 2001| . 

To precisely define the knee dislocations and disclinations, caused by the S'£'(3)— jolt (F,T), we 
first define the coordinate co- frame, i.e., the set of basis 1-forms {da;'}, given in local coordinates 
x^ = (x^ , x^, x^) = {x, y, z), attached to the moving segment's center-of-mass. Then, in the coordinate 
co-frame {dx^} we introduce the following set of the knee plastic-deformation-related S'-E(3)— based 
differential p— forms (see [Ivancevic and Ivancevic 2006c[ ITvancevic and Ivancevic 2007e] ): 

the dislocation current 1-form, J — Jidx^; 

the dislocation density 2-form, a — 2'^ij dx^ A dx^ ; 

the disclination current 2-form, S = o'^y '^^* ^ dx^ ; and 

the disclination density 3~form, Q = -^Qijk dx^ A dx^ A dx^ , 
where A denotes the exterior wedge-product. According to Edelen [Edelen 1980[|Kadic and Edelen 1983| . 
these four S* £^(3)— based differential forms satisfy the following set of continuity equations: 

a = -dJ - S, (6) 

Q = -dS, (7) 

da = Q, (8) 

dQ = 0, (9) 

where d denotes the exterior derivative. 

In components, the simplest, fourth equation ^, representing the Bianchi identity, can be rewrit- 
ten as 

dQ = diQ[ijk] dx^ A dx^ A dx^ A dx'' = 0, 

where di = d/dx'^, while ^[y...] denotes the skew-symmetric part of %.... 
Similarly, the third equation ([5]) in components reads 

TT.Qijk dx^ A dx-' Adx~ — dkOi[ij] dx A dx^ A dx^ , or 

Qijk = -69fea[y]. 
The second equation ([7]) in components reads 

—Qijkdx' Adx^ Adx'' = -dkS[ij]dx'' Adx' A dx\ or 
Qijk — 69fe5[ij]. 
Finally, the first equation ^ in components reads 

-dij dx^ A dx-' = {djJi ~ -Sij) dx^ A dx-' , or 



— ^^j^i ^ij • 



In words, we have: 



• The 2-forni equation ([5]) defines the time derivative d —-^ctij dx^ f\dx^ of the dislocation density 
a as the (negative) sum of the disclination current S and the curl of the dislocation current J. 
This time derivative is caused by the translational part of the SE(3) jolt. 

• The 3-form equation ([7]) states that the time derivative Q =-^Qijk dx^ A dx^ A dx^ of the discli- 
nation density Q is the (negative) divergence of the disclination current S. This time derivative 
is caused by the rotational part of the SE(3) jolt. 

• The 3-form equation ^ defines the disclination density Q as the divergence of the dislocation 
density a, that is, Q is the exact 3-form. 

• The Bianchi identity ^ follows from equation ([5]) by Poincare lemma [Ivancevic and Ivancevic 2006c[ 
llvancevic and Ivancevic 2007ej and states that the disclination density Q is conserved quantity, 
that is, Q is the closed 3-form. Also, every 4-form in 3D space is zero. 

From these equations, we can conclude that the knee dislocations and disclinations are mutually 
coupled by the underlaying S'i?( 3)— group, which means that we cannot separately analyze transla- 
tional and rotational knee injuries. This result supports the validity of the combined loading hypoth- 



3 Conclusion 

Based on the previously developed covariant force law |Ivancevic and Ivancevic 2006a[|Ivancevic and Ivancevic 2007ej . 

its recent application to traumatic brain injury [Ivancevic 2008] . and using as an example a human 

knee joint, in this paper we have formulated a new coupled loading-rate hypothesis for the generic 

musculo-skeletal injury. This new injury hypothesis states that generic cause of human joint injuries 

is an external S'i?(3)— jolt, an impulsive loading hitting a joint in several degrees-of-freedom, both 

rotational and translational, combined and simultaneously. To demonstrate this, we have developed 

the vector Newton-Euler mechanics on the Euclidean iS'£'(3)— group of the knee micro- motions. In 

this way, we have precisely defined the concept of the 5'i?(3)— jolt, which is a cause of two kinds of 

rapid joint discontinuous deformations: (i) mild rotational disclinations, and (ii) severe translational 

dislocations. Based on the presented model, we argue that we cannot separately analyze localized 

joint rotations from translations, as they are in reality coupled. To prevent human musculo-skeletal 

injuries we need to develop the m,usculo- skeletal SE(3)-jolt awareness, e.g. never overload a fiexed 

knee, avoid any kind of collisions. 
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